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We study the ground state of a d-dimensional Ising model with both long range (dipole-like) 
and nearest neighbor ferromagnetic (FM) interactions. The long range interaction is equal to r~^ , 
C/3 ' p> d, while the FM interaction has strength J. If p > d + 1 and J is large enough the ground state 

is FM, while ifd<p<d+l the FM state is not the ground state for any choice of J. In d = 1 we 
Cy ■ show that for any p > 1 the ground state has a series of transitions from an antiferromagnetic state 

S ' of period 2 to 2/i-periodic states of blocks of sizes h with alternating sign, the size h growing when 

I , the FM interaction strength J is increased (a generalization of this result to the case < p < 1 is 

' O ' also discussed). In d > 2 we prove, for d < p < d+ 1, that the dominant asymptotic behavior of the 

C ' ground state energy agrees for large J with that obtained from a periodic striped state conjectured 

O , to be the true ground state. The geometry of contours in the ground state is discussed. 

Cq ; I. INTRODUCTION 
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, ^ . There has been and continues to be much interest in niesoscopic pattern formation induced 

vQ ' by competing short and long range interactions [1-4]. Of particular interest is the competition 

tH- . between dipole-like interactions bewteen spins which decay as r^'''+^\ d the dimension of 

(^ ' the lattice, and short range exchange interactions [5, 6]. For d = 2, this type of competitive 

^^D , interaction is believed responsible for many of the observed patterns in thin magnetic films [7] 

^P ■ and other quasi-two dimensional systems, including Langmuir monolayers [8], lipid monolayers 

[9], liquid crystals [10], polymer films [11] and two-dimensional electron gases [12, 13]. 

The simplest models to describe such systems are Ising spins with a nearest neighbor fer- 
romagnetic interaction and a power law long range antiferromagnetic pair potential, of dipole 
(or Coulomb) type [6, 14-19]. The zero temperature phase diagram of these models has been 
pH ' thoroughly investigated over the last decade and a sequence of transitions from an antiferromag- 

Q , netic Neel state to periodic striped or lamellar phases with domains of increasing sizes has been 

O ' predicted, as the strength of the ferromagnetic coupling is increased from zero to large positive 

values. These theoretical predictions are mostly based on a combination of variational tech- 
niques and stability analysis: they start by assuming a periodic structure, proceed by computing 
the corresponding energy and finally by comparing that energy to the energy of other candi- 
date structures, usually by a combination of analytical and numerical tools. These calculations 
give an excellent account of the observed "universal" patterns displayed by the aforementioned 
quasi-two dimensional systems. However they run the risk of overlooking complex microphases 
that have not been previously identified [20]. This risk is particularly significant in cases, as 
those under analysis, where dynamically (e.g. in Monte Carlo simulations) the domain walls 
separating different microphases appear to be very long lived as the temperature is lowered [14]. 
In order to settle the question of spontaneous pattern formations on more solid grounds it 
would be desirable to be able to first prove periodicity of the ground state and then proceed 
with a variational computation within the given ansatz. The problem is not simple. Most 
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of the mathematically rigorous techniques developed for obtaining the low temperature phase 
diagram of spin systems, e.g. the Pirogov-Sinai theory [21], depend on the interaction being 
short range. Only methods based on reflection positivity [22] or on convexity [23-27] seem 
applicable to the kind of systems considered here. In this paper we apply reflection positivity 
to give a characterization of the ground states of the one-dimensional system and to give in 
this case a full justification of the variational calculation based on the periodicity assumption. 
Moreover we obtain rigorous upper and lower bounds on the ground state energy in higher 
dimensions, thus providing a quantitative estimate of the possible metastability of the striped 
or lamellar phases. 



A. The model 

Given an integer N, let An be a "simple cubic" d-dimensional torus of side 2N and a G 
|_l_j^|Ajv_ Any configuration ct is a sequence of (7^ = ±1 labelled by i e Ajy. The Hamiltonian 
of this Ising model with periodic boundary conditions is taken to be 

d ^ 

Hn{2.) = -JY.Y1 ^^^«+-^ +11'^^ MJ - i) <^j . JpU - «) = H I ■ _ ■ ^ 2nN\P ^^'^^ 

iGAjv fe=l (i,j) nez'^ ' ' 

where in the first sum e^ is the unit vector in the fc-th coordinate direction (and if z + e^; ^ An 
we define o'i-i-e^ = <Ji-2Neh)i the second sum runs over generic distinct pairs of sites in An 
and p > d. The first term in (1.1) will be called the exchange energy term and we will take 
J > while the second one will be called the dipolar energy term. In d > 1 we shall define 
|z — j| to be the usual Euclidean distance between i and j. Note that the sum over n in the 
definition of Jp{j — i) makes sense as long as p > rf. However some of the results discussed 
below hold true even in the case d—l<p<d, provided the definition of Jp(j — i) is replaced by 
Jp{j -i) = \j - i\-P + (27V - \j - z|)-P, VI < \j -i\<2N ~-\. In the following we wiU restrict 
attention to the case p > d and we will comment on possible generalizations of our results to 
d~ 1 < p < din the Remark after Theorem 2. See also the Remark after Theorem 3. 



B. Main results 

We wish to determine the ground state of model (1.1) for different values of p and J. A 
first remark is that in any dimension for p > d -\- 1 and J large enough the ground state is 
ferromagnetic and is unique modulo a global spin flip. This is, in d > 1, a corollary of the 
contour estimates in [28], and we shall reproduce its proof as a byproduct of our analysis. For 
this case the low temperature Gibbs states are also known: for d—1 and p > 2 there is a unique 
Gibbs state for any (3 < +00 [29] while for d>2, p>d+l and /3 large enough, there are two 
different pure states obtained as the thermodynamic limits of the equilibrium states with + or 
— boundary conditions [28]. On the contrary, for d < p < d + 1, the ferromagnetic state is not 
the ground state in any dimension, for any value of J. Instability of the ferromagnetic state 
in presence of a long range antiferromagnetic interaction of the form 1/r^, d < p < d + I, was 
first noted in [30] . The reason for this instability lies in the divergence of the first moment of 
the long range interaction which makes the ferromagnetic state unstable towards flipping spins 
in a large enough domain. 

Our main results give a characterization of the ground state in the whole regime p > d. 
For d = 1 the characterization is complete, in the sense that for any p > 1 we can compute 
the ground state energy per site and we can prove that the corresponding ground state 
configurations are periodic and consist of blocks of equal size and alternating signs (a block is 
a maximal sequence of adjacent spins of the same sign). 

Theorem 1. (Ground state energy, d=l). 

For d = 1 and p > I, the ground state energy Eo{N) of the Hamiltonian (1.1) satisfies 

lim -±iJo(iV)=mine(/i) (1.2) 

N^oo Zl\ /ifEZ+ 



where e{h) ~ liniTv^oo (2 A^ /i) ^Eper{Nh) and Eper{Nh) is the energy of a periodic configura- 
tion on a ring of 2Nh sites consisting of blocks of size h with alternating signs. 

The function e{h) is explicitly computable, see Section III. The minimization of e(/i) with 
respect to h can be performed exactly and, as expected, for p > 2 and J larger than an explicit 
constant Jq = T{p)^^ L daaP^^e^'^{l — e^")^^, the minimizer is ft- = +cxd (i.e. the ground 
state is ferromagnetic). For 1 < p < 2 and all J > or p > 2 and < J < Jq the minimum 
of e{h) over the positive integers is attained at an integer h*{J) which is a piecewise constant 
monotone increasing function of J. Thus there is a unique integer minimizer of e(/i) for almost 
all values of J which jumps at a discrete set of values of J. At those values of J the minimum 
is attained at two consecutive integers h*{J) and h*{J) + 1 and the corresponding ground 
states, for suitable values of N, are the configurations with all blocks of the same size, either 
h*{J) or h*{J) + l. 

Theorem 2. (Finite volume ground states). 

Let d— 1 and p > I. At the values of J such that the minimum of e{h) is attained at a single 
integer h*{J), then in a ring of length 2N , such that N is divisible by h*{J), the only ground 
states are the periodic configurations consisting of blocks of size h*{J). At the values of J such 
that the minimum is attained at two consecutive integers h*{J) and h*(J) + 1, then in a ring 
of length 2N , such that N is divisible both by h*{J) and h*{J) + 1, the only ground states are 
the periodic configurations consisting of blocks of the same size, either h*{J) or h*{J) + 1. In 
both cases there is a finite (independent of N ) gap between the energies of the ground states 
and of any other state. 

Remarks. 

1) Stability for < p < 1. The proofs of Theorems 1 and 2 do not really require p > \ and, 
provided we change the definition of Jp(j — i) as described after (1.1), they can be easily 
adapted to cover the cases < p < 1. In particular for any p > Q the ground state energy 
can be computed by minimizing energy among the periodic states of blocks with alternating 
signs and the computation shows that the ground state energy is extensive (i.e. it scales 
proportionally to iV as iV ^ oo), even for < p < 1 (a case in which thermodynamic stability 
is not a priori guaranteed). 

2) Infinite volume ground states. A corollary of Theorem 2 is that the infinite volume periodic 
configurations g_j G {±1}^ with blocks all of the same size h*{J) (or h*{J) + 1, if J is a value at 
which the minimum of e(/i) is not unique) and alternating sign arc infinite volume ground state 
configurations, in the sense that they are stable against bounded perturbations [21]; i.e. given 
a finite set X C Z, if Rx is the operator hipping the spins in X, the (finite) energy difference 
between Rx2.j and CTj is positive. This can be proven by writing the energy difference between 
the infinite configurations Rx2Lj and g_j as the limit of the energy differences between the finite 
volume configurations on rings of length 2N (N divisible by h*{J)) obtained by restricting 
RxQ_j and g_j to the finite rings. By Theorem 2, the differences between the energies of the 
finite volume approximations of RxSLj and CTj are positive. This implies that CTj is an infinite 
volume ground state configuration. 

A stronger result which follows from our analysis is that for any J > we can choose a 
sequence of integers {TVijign, Ni ^^^> oo, such that the ground states on the rings of lenghts 
2Ni are periodic configurations of blocks of size h*{J) and the corresponding sequence of 
ground state energies per site eo{Ni) is the sequence of "lowest possible specific energies", i.e. 
any other sequence {A^/jigN, N^ -^^ — > oo, has eo{Nl) > eo{Ni), eventually as i ^ oo. 

In dimension larger then one we do not find an exact asymptotic expression for the 
ground state energy. Still, for large J, wc find rigorous upper and lower bounds to the ground 
state energy. The result is the following. 

Theorem 3. Let d > 2. If p = d + 1 and J is large enough, there exist positive, 
J -independent, constants Ci,C'2 such that 

epMi-J) - Cie-^l=^^-il"' < lim -^Eo{N) < epMi-I) - Cae-^l^^'-il"' , (1.3) 



where epMiJ) is the energy per site of the ferromagnetic configuration {a)i = +1 and ISd-il is 
the volume of the d — 1 dimensional unit sphere. Ifd<p<d+l and J is large enough, then 
there exist positive, J -independent, constants Ki,K2 such that 

epMiJ) - KiJ^^^S^ < lim — — -Ea{N) < epuiJ) - K2J^^^^^ . (1.4) 
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Remark. {The case p < d). Contrary to the proofs of Theorem 1 and 2, the proof of Theorem 
3 crucially relies on the summability of the potential, i.e. on the condition d < p < d + I. It 
would be of great interest to extend the proof to the case d — 1 < p < d and to the case p = I 
in dimension d = 2,3 (Coulomb case). A thorough discussion of the case d— l<|7<(i+lin 
dimension d > 2, including results related to those of Theorem 3, is provided by Spivak and 
Kivclson [13]. The special case p = 1 with d = 2 was treated in [16], and in full generality in 
[31]. The Coulomb case p = 1 in d = 3 was considered in [17]. 

The upper bounds in (1.3) and (1.4) follow from a variational computation. The best 
known constants C2 , K^ are obtained by minimizing over the periodic striped configurations 

{^strinedii ~ (—1)'*^^^', where ii is the first component of i and [x] is the largest integer 
less than or equal to x. It is remarkable that, for J -^ cxd, the minimum over the striped 
configurations provides a variational energy that is lower than the one obtained by minimizing 
over the periodic checkerboard configurations; this was shown in [14] for d — 2 and p — 2>. K 
computation analogous to that in [14] allows to prove a similar result in higher dimensions 
for d < p < d + 1. This result together with numerical studies support the conjecture 
that the ground state of (1.1) in d > 2 and d < p < d -\- 1 is m fact a periodic striped 
configuration (at least asymptotically for large J). In this paper we prove the lower bounds in 
(1.3) and (1.4), which show that for large J the difference between the ground state energy 
and the ferromagnetic energy scales exactly as predicted by the above mentioned variational 
computation. The proof of the lower bound is based on an energy argument which also gives 
estimates on the allowed shapes and sizes of Peierls' contours in the ground state. 



II. OUTLINE OF THE PROOFS 

Before presenting the proofs of Theorems 1, 2 and 3 in detail, we give a short outline. 

The proofs of Theorem 1 and 2 are based on reflection positivity. It is known [32] that for 
J = the Hamiltonian (1.1) is reflection positive in any dimension. In d = 1 this means that, 
given any configuration of spins on a ring of 2A^ sites g_ ~ (ct^at+i, . . . , cr_i, ctq, fxi, (J2, ... , cn) 
and a bond b = (i,i + 1) connecting site i with its neighbor i + 1, the average energy 
of the two conflgurations obtained by reflecting around b is always lower than the original 
one; e.g. choosing b — (0, 1) the average of the energies of the two reflected configura- 
tions (-(Tat . . . , -0-2, -O-i, (Ti, (72, . . . , CTAf) and (cT-AT+i . . . , cr-i, cro, -cro, -0--1, . . . , -CT-Af+l) is 

always lower than the energy of the original configuration a_. Using this property and repeat- 
edly reflecting around different bonds, one can prove that for J = the ground state is the 
antiferromagnetic alternating state. Note that periodic boundary conditions are necessary for 
repeated reflections around different bonds. 

Unfortunately, as soon as J > 0, the Hamiltonian (1.1) is not reflection positive anymore. 
Still, one can think of making use of the reflection symmetry around bonds b = {i,i + 1) 
separating a spin ai from a spin CTi+i = — (7^: such reflection does not change the exchange 
energy between Oi and (Ji+i and lowers the dipole energy. By repeatedly reflecting around such 
sites one could expect to be able to reduce the search for the ground states just to the class 
of periodic configurations of blocks of the same size and alternating sign and explicitly look 
for periodic configuration with minimal energy (we recall that by block we mean a maximal 
sequence of consecutive spins all of the same sign). However there is a difficulty: because of 
periodic boundary conditions, in order not to increase the exchange energy in the reflection, 
one should reflect around bonds b not only separating a -I- from a — spin, but with the further 
property that the bond b' at a distance N from b also separates a + from a — spin: but in a 
generic configuration a there will be no bond b with such a property! 



A possible solution to this problem is to cut the ring into two uneven parts both containing 
the same number of blocks but not necessarily of the same length. The configurations obtained 
by reflecting the two uneven parts will have a lower energy but in general different lengths 
2N',2N". The idea is to reflect repeatedly in this fashion, keeping track of the errors due 
to the fact that the length of the ring is changing at each reflection. A convenient way of 
doing this, exploited in Section Iff, is to rewrite the spin Hamiltonian (1.1) as an effective 
Hamiltonian for new "atoms" of "charges" hi , corresponding to the spin blocks of size hi . The 
new effective Hamiltonian E(. . . , /i_i, ho, hi, h2, . . .) is again reflection symmetric (in a slightly 
different sense, though) and its explicit form allows for an easy control of the finite size errors 
one is left with after repeated reflections. Some technical aspects of the proofs of Theorem 1 
and 2 are given in the Appendices. 

The proof of Theorem 3, in particular the lower bounds in (1.3) and (1-4), is based on a 
Peierls' contour estimate described in Section IV. For large ferromagnetic coupling J, we shall 
describe the ground state configuration(s) in terms of droplets of — spins surrounded by + spins 
and of contours separating the + from the — phases. The requirement that the energy of the 
ground state configuration is minimal imposes bounds on the geometry and the energy of such 
droplets, implying the lower bound in (1.4). Contrary to the methods exploited in the proofs of 
Theorem 1 and 2, the proof of Theorem 3 is robust under modifications both of the boundary 
conditions and of the specific form of the interaction potential. 



III. ONE DIMENSION 

In this section we want to prove Theorems 1 and 2. Restricting to d = 1, we shall consider 
any configuration g_ S {±1}^™ on the ring of length 2N as a sequence of blocks of alternating 
sign, where by definition a block is a maximal sequence of adjacent spins of the same sign. We 
note that, due to the periodic boundary condition, the number M of blocks on the ring is either 
1 (if the state is ferromagnetic) or an even number. We shall denote hy hi, i = \, ... , M , the 
sizes of such blocks. A block will be called a + block (— block), if its spins are all of sign +1 

(-!)■ 

We want to prove that the sizes hi of the blocks in the ground state are all equal, at least 
for N large enough. The strategy will be the following. Given any configuration ct^ in the 
ring Km with M > 2 blocks of alternating signs of sizes h_M^i, ■ ■ ■ ,hA£, we will rewrite the 
energy HN{q_j^) as a function of M and the ft-i's, i.e. -ffAr(CT^) = EN{M,{h]^,hj^)'j, where 
hj^ — {h_M__^_■^^, . . . ,ho) and hp^ = {hi, . . . , /im ). Setting h — {hL,hfi), we will show that, for M 
and N fixed, the Hamiltonian Em{M,Ii) is reflection positive with respect to the reflection 

^, , MM , ^ 

Oh^ = /i_j+i , — — <i < — , (3.5) 



that is 



EN{M,h) > i [EN{M,{h^,0h^)) + EN{M,{9h^,h^))] , (3.6) 



where 9h]^ = {ho, . . . , h_M+i) and 0hj^ — {hM_, . . . , hi), with of course the signs of the spins 
in the reflected blocks being opposite to what they were originally. Using repeatedly this 
symmetry, we will get bounds from above and below for the ground state energy, in terms of 
the energy of periodic configurations. 



A. The integral representation 

In order to show reflection positivity of EN{M,h), e.g. (3.6), we need to look more closely 
at the structure of the Hamiltonian. A straightforward calculation in Appendix B gives the 
energy Ej^{M,h) of a configuration in An of M blocks of sizes h_ivL_^_^, . . . , hM_ (recall that M 
is even and p.b.c. are assumed) as: 

/■°° da ■="" 

EN{M,h) = -2J{N - M) + 2N I —— a^-^ 



T{p) (l-e-")2(l 



•{ J2 [/i*(l-e-")(l-e-2"W)-(l-e-"'*0(l-2e-"(2^-"')+e-2"^^)] + 






-ah. 



)(i-e-"''o[ n ^ 

i<k<j 



ahk 



Yl £-"'''=]}. (3.7) 



-■M.<k<i 



The expression in braces in the last two hnes of (3.7) can be pictoriaUy interpreted as the 
energy of a system of particles on a ring of M sites with hi>\ particles at site i. There is an 
on site energy as in the second line of (3.7) and a "many body" energy as in the third line. 

The energy in (3.7) has the remarkable property of being reflection positive with respect to 
the reflection (3.5). This can be proven by rewriting the expression in braces in the last two 
lines of (3.7) as 



/(a) \HR{a, hji) + Hnia, eh^^) - ^ F,{a, h^) F^{a, ^/t^) 



da 



(3.8) 



where 



v{a) = 



vP-i. 



r(p) (l-e-")2(l-e-2"A^) ' 

M 
2 



-ahk 



-ahk 



F^{a,h^) = J2{-iy{l-e-^'^') Yl 

i=l l<fc<i 

M 

i^2(a,/.^) = ^(-1)^(1 -e-"'^0 n 
»=i i<fc<f- 

Let us define H(a,h) as follows: 

2 

H{a,h) = HR{a,h^) + Hnia, §h^) - ^ F,{a,hj^) F,{a, Oh^) 



(3.9) 



(3.10) 



so that EN{M,h) = -2J{N - M) + 2N ^^ ^ aP-^ ^^^-^ + /g°° dav{a)H{a,h). 

Using the fact that Hji{a,hj^) — F[ji{a,6hji) and that, by Schwarz inequality, 
F,{a,hji)F,{a,§h^) < l[F,{a,hi^)F,{a,hji) + F.iajh^) F^{a,eh^% we find 



da via) H(a, h) = 

2 

dav{a) Hu^a.hji) -\- HR{a,Oh^) - y^ f ^ (a, /ij^) f^ (a, Oh^ 



(3.11) 



> 



> - / dai^ia) [H{a,{hL,ehL)) + H{a,{ehji,h^)) 

so that (3.6) is proved. 

Rcficcting repeatedly with respect to different bonds we end up with 

M 

/dai'{a) H{a,h) > — > dai'{a) H(a, {hi,hi, . . . ,hi,hi)) 



(3.12) 



and EN{M,h) > ttX)- Mj_i ^N{M,{hi,hi, . . . ,hi,hi)) , which is an example of the chess- 
board inequality, see Theorem 4.1 in [22]. 

B. Bounds for 1 < p < 2 

We now temporarily restrict to the case 1 < p < 2. A key remark is that in this case there 
exists a p-dependent constant Kp, explicitly computable as described in Appendix A, such that 
the blocks in the ground state have sizes hi all satisfying the following apriori bound: 

hi < 12e-' , if p = 2 /"3 i o^ 

/i, < i^p ji/(2-p) if l<p<2, ^'^•^•^^ 

Eq. (3.13) shows in particular that, if 1 < p < 2, the FM state is not the ground state for any 
finite J when N is sufficiently large; in particular in the ground state the number M of blocks 
is M > 2 (and necessarily even). Note also that since all blocks satisfy 1 < hi < h^ax, with 
hmax given by (3.13), it must be M < 2N < h^axM, that is N and M are of the same order 
as iV ^ oo. 

To compute the ground state energy, let us introduce the auxiliary partition function 



Qj^^Y. e-^^"(^^^) (3.14) 



where the * means that we are summing only over the choices of M and h compatible with 

M 

the constraint X) -_m , , hi ~ 2N and with the bounds (3.13). We can then obtain the ground 

state energy per site eo(A^) by taking the hmit — hm^g^oo -h logQAr — 2Neo{N). 
Using (3.12) we get: 

M 

Qn<J2 n e-^^^"(^^'('^--''')). (3.15) 

M,h i=-M + i 

By (3.7), wc find that ■^En{M, {h, . . . , h)) can be written explicitly as 

-E^{M,ih,...,h)) = -Jh + Ah + 2J + j^ ______|_2tanh— + 

-ah\2 



1 (1 - e-"'') 



_L(2e-"(2N-/.) ^ ^-2»N _ g-ah(M-l)^| (3^^g) 



(1 - e-2"W) I ^ p- 
where the constant A in the r.h.s. is A — [^ t^ aP^^jr-^ — -rj. Note that the absolute value 

Jo r(p) (l-e-°) 

of the term in the last line can be bounded above by KN^^, for a suitable constant K. Then, 
recalling the definition of EpJr{Nh) (see Theorem 1) and defining e{h) as 

WC find that the r.h.s. of (3.15) and thus Qn can be bounded above by 

M 

Qn<J2 n e-^'^'^'^^^ef^'"''''' (3.18) 

Using the fact that the number of terms in the sum is less than 2^^, we can bound the r.h.s. of 
(3.18) from above by 2'^'^e^^'^^''''' e''^'^^^" , where cq = min;,ez+ e{h), so that we find eo(A^) > 
eo - K'/NP. 



It follows from the explicit expression of e{h) given in (3.17) that niin;jg2+ e{h) — e(h*(J)), 
where the integer h*{J) is a piecewise constant monotone increasing function of J . Thus there 
is a unique integer minimizer of e{h) for almost all values of J which jumps at a discrete set 
of values of J. At those values of J the minimum is attained at two consecutive integers h*(J) 
and h*{J) + 1. As discussed below, the corresponding ground states are, asymptotically for N 
large, the configurations with all blocks of the same size, either h*{J) or h*{J) + 1. If ft, is not a 
minimizer, then e{h) is separated by a gap from eo- These properties follow from the fact that 
dhe{h) — has a unique solution for /i e M and the solution is a local minimum. 

Let us now turn to the problem of finding a lower bound for Qn- If N is divisible by h* , 
then of course 

SO that we find 

eo{N)^eo + 0{^). (3.20) 

If A'' is not divisible by h* the error in the r.h.s. is replaced by 0(N^^). This follows from the 
fact that we can bound Qn from below by restricting the sum in (3.14) to a configuration with 
all but one block of sizes h*: the last block being of size h* < h < 2h*. The energy of such 
configuration is 2Ne(h*) + 0(1). This concludes the proof of Theorem 1 in the case p < 2. 

C. The case p > 2 

Let us now discuss the case p > 2. It is straightforward to verify that in this case, if 
J > Jo = /g j4^a^~"'" /-, _^.-a-):j , then the minimizer of e(/i) is h = +cxd. Correspondingly we 
can prove that the ground state of (1.1) is the ferromagnetic state. In fact let us assume by 
contradiction that in the ground state there is a block Bh of finite size h. Then it must be 
true that the energy AE{h) needed to reverse the sign of all the spins in _B^ is nonnegative: 
AE{h) > 0. On the other side AE{h) < -4J + 2Ei{h), where Ei{h) is the dipolc energy 
between the + block Bh and an external sea of + spins. It is straightforward to check that 

and this leads to a contradiction. This proves Theorem 1 for p > 2 and J > Jo- 

If on the contrary J < Jq, a repetition of the proof in Appendix A leads to the bound 
hi < K2/{Jo — J)^/(P^^) on the sizes hi of the blocks in the ground state. Then we can repeat 
the proof above to get the desired result (1.2) and this concludes the proof of Theorem 1. ■ 



D. Uniqueness of the ground state 

In this subsection we prove Theorem 2. If p > 2 and J > Jq the statement is a corollary 
of the proof above: in fact the contradiction obtained after (3.21) shows that in this case the 
ferromagnetic state is the unique ground state. 

Let us then consider the cases l<p<2orp>2 and J < Jq, in which wc have an apriori 
upper bound on the sizes of the blocks. Let us first consider the case in which the minimizer 
h* of e{h) is unique and let us choose 2A^ = Mh* . Let ft be a configuration in An for which 
the set X = {i : ft" 7^ ft*} is non empty. Note that the number of blocks M° in h is not 
necessarily equal to M. Using again the chessboard estimate (see Theorem 4.1 of [22]), the 
energy EN{M^,h^) of the configuration ft° can be bounded below as: 

EN{M',h')>^ V £;^(M°,(ftO,...,ft^))> V (ftOe(ft°)-^) (3.22) 



-+i 



Using (3.20) we find that 

EN{M\lf) 2Neo{N) > ^ h'^{e{h'^) e{h*)) + 0(-^) , (3.23) 



iei 



so that, since e(/i") — e{h*) > Ae, we have that for A^ big enough the unique ground state for 
2N — Mh* is the state with M blocks all of sizes h* . Similarly, if 2N is not divisible by h* , 
any configuration with a sufficiently big number of "wrong" blocks (i.e. of blocks with sizes 
^ h*), will be separated by a gap from the ground state. 

Finally, let us consider the case in which there are two minimizers h*,h* + 1. The 
same proof as above goes on to show that, if N is divisible either by h* or by /i* + 1, then 
the ground state energy is given by (3.20) and the energy of configurations with blocks of 
"wrong" sizes (i.e. different both from h* and h* + 1) is asymptotically larger than the energy 
of the ground state. Similarly, if N is not divisible either by h* or by ft.* + 1, the energy of 
configurations with a sufficiently big number of blocks with wrong sizes will be larger than the 
ground state energy. 

We are left with considering the case of a configuration with blocks' sizes all equal ei- 
ther to h* or to ft* + 1. Let us denote by g_{h*) the 2ft*-periodic configuration with all blocks 
of size ft* and by g_{h* + 1) the 2(ft* + l)-periodic configuration with all blocks of size ft*. Can 
a configuration with a finite fraction of blocks of sizes ft* and ft* + 1 be a ground state? 

Let us consider a configuration ft with ft? G {ft*, ft* + 1} and the total number of blocks equal 
to M in a volume 2N , with N divisible by ft* or ft* + 1 and M divisible by 4 (this condition on 
M is not restrictive: if M were not divisible by 4 we could consider the configuration obtained 
by doubling the system). Let 

I, "H {» : /,o ^ /,o^^ =./,*}, I, "H {^■.h'^^ hi, = ft* + 1} , J2 "^^ {z : ft? 7^ ft°+i} , (3.24) 

where of course ftM. , , = h^_M_,,- Note that \Iq\ + |22|/2 = Mo, where Mq is the number of 

blocks of size ft*, and \Ii\ + |22|/2 = Mi, where Mi is the number of blocks of size ft* + 1. 
Using again the chessboard inequality, we find that 

M 

EM{M,lf)>l^ E i^A^(M,(ftO,ft°+i,ft°+i,ft°,...,ft°,ft°+i,ft°+i,ft°)) (3.25) 

Now, except for an error of 0(M^(p^^'), the r.h.s. of (3.25) can be rewritten as |Xo|ft*eo + 
|Xi|(ft* + l)eo + 1X2 1 (ft* + \)e, where e is the (infinite volume) specific energy of the (4ft* + 2)- 
periodic configuration obtained by repeating periodically the configuration (ft*, ft* + 1, ft* + 1, ft*) 
over the volume. The key remark is that e — cq — 5 > Q: this is proven in Appendix C. Then, 
using the fact that |Xo| + |22|/2 = Mq (with Mq the number of blocks of size ft*) and that 
\1i\ + 1^2 1/2 — Ml (with Ml the number of blocks of size ft* + 1), we get 

EN{M,hl-2Neo{N)> |J2|(ft* + ^)(e - eo) (3.26) 



2' 



and the proof of Theorem 2 is concluded. 



IV. HIGHER DIMENSIONS 

In this section we want to prove Theorem 3. As already remarked after Theorem 3, for 
d — 2 and p — i the upper bound was obtained in [14] by minimizing over the periodic striped 
configurations (and checking that asymptotically for large J such an upper bound is better 
than the one obtained by minimizing over the periodic checkerboard configurations). The 
upper bound in dimension higher than 2 and for d < p < d + 1 can be again obtained by a 
variational computation, minimizing the energy over the periodic striped configurations. We 
do not repeat the details here; the result is provided by the upper bounds in (1.3) and (1.4). 
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We now focus on the lower bounds in (1.3) and (1.4). We need to introduce some definitions; in 
particular via the basic Peierls construction we introduce the definitions of contours and droplets. 
Given any possible configuration a on An, and in particular the ground state configuration(s), 
we define A to be the set of sites at which (Ti = —1, i.e. A = {zsA : ai = —1}. We draw 
around each i G A the 2d sides of the unit (rf — 1) dimensional cube centered at i and suppress 
the faces which occur twice: we obtain in this way a closed polyhedron r(A) which can be 
thought as the boundary of A. Each face of r(A) separates a point « £ A from a point j ^ A. 
Along a. d~2 dimensional edge of r(A) there can be either 2 or 4 faces meeting. In the case of 
4 faces, we deform slightly the polyhedron, "chopping off" the edge from the cubes containing 
a — spin. When this is done r(A) splits into disconnected polyhedra 71, ... ,7^ which we shall 
call contours. Note that, because of the choice of periodic boundary conditions, all contours 
are closed but can possibly wind around the torus Ajy. The definition of contours naturally 
induces a notion of connectedness for the spins in A: given i, j £ A we shall say that i and j are 
connected iff there exists a sequence (i — iQ,ii, . . . ,in — j) such that im, im+i, m = 0, . . . , n~ 1, 
are nearest neighbors and none of the bonds (i^, im+i) crosses r(A). The maximal connected 
components 5i of A will be called droplets and the set of droplets of A will be denoted by 
I?(A) = {Si, . . . , Ss}. Note that the boundaries r{Si) of the droplets di £ 2?(A) are all distinct 
subsets of r(A) with the property: Lll^-^r{Si) = T{A). 

Given the definitions above, let us rewrite the energy in (1.1) as 

HNi^) = \AN\eFMiJ) + 2J J2 1^1- J2 Edrp{5) , (4.27) 

7Gr(A) <5eX'(A) 

where epMi-J) is the energy per site of the ferromagnetic configuration ai = +1 and 

de f 

Edip{5) — 2 X^iGcj '^ieA" '-^pi^ ^ .?)• ^^^ ^^ arbitarily choose a number £ > 1 (to be conveniently 
fixed below) and let us correspondingly rewrite Edip{6) as 

Ed,p{S)^2j2J2 n\^-J\<e)M^-J) + 2J2J2 ]l(K-j|>^)Jp(*-j) (4.28) 

ies jeA" ies jeA" 

Denoting the last term by E^^{5), we note that X^aeprA) -^^^('^) '^^'^ ^^ bounded above by 
2min{|A|,|A=|}$p(^) < |A|$p(£), where $p(£) = E|„[>£ 1"!"^ < const.£-^P-'^\ For large £ 
the constant is smaller than 2d\^d\/ {p — d) where \Sd\ is the volume of the d-dimensional unit 
sphere. The first term in the r.h.s. of (4.28) can be rewritten as 

i?^^(<5) = 2^ Jp(n)^^]l(^-J.= n)<2^^^ J2 M^ - J = n) (4.29) 
|n|<£ ieS jeA" \n\<e'^' leS jeid\s 

where in the last inequality we neglected an error term vanishing in the thermodynamic limit. 
Now, the number of ways in which n ~ (rii, . . . , Ud) may occur as the difference i — j or 

j — i with i e (5 and j ^ i5 is at most X^-Yerd?) Xi=i iTM'^il; where |7|i is the number of faces 
in 7 orthogonal to the z-th coordinate direction. To prove this, draw a path on the lattice of 
length |ni| + • • • + |nrf| connecting i and j. Such a path must cross a face of 7 and if this face 
is orthogonal to the z-th coordinate axis, can do so only in \ni\ ways and the claim follows. 
The conclusion is that 



-ieT(S)\n\<l' '' i=l 76r(,5) \n\<l 



„|P /l^l MM -M Z_^ I M ,^^ |^|p 



The sum X|n|<£ l^'il/l"!'' '^^^ ^^ bounded above by j^^,.^^d'^x\xi\/\x\P + const, where the 
constant is smaller than 2pd\Eid\2'P~'^ / {p — d) and jS^d is the volume of the d-dimensional unit 
sphere. Then (4.30) implies E'^^{6) < 2-^p{l) + 2pd\Sd\2P-'^ /{p- d), where *p(f) = |Sd-i|fog^ 
iip = d+l and*p(^) = \^d-i\{(-'^^^^^ ~ l)/ {d + I - p) iid<p<d+l. 
Putting these bounds back into (4.27), we find: 

E»(«) > |A„|e„, + 2 (. - ?f^ - *,(.,) i: H - a;5MM , ,,3:, 

^ ^ / 7er(A) ^ 



11 



where A is the set of — sites in the (unknown) ground state configuration. Choosing i in such 
a way that *p(£) = J - pd\Sd\2P-'^ / {p - d), wc find the lower bounds in (1.3) and (1.4). This 
concludes the proof of Theorem 3. ■ 



As a side remark, let us note that a discussion similar to the one above implies that iip > d+1 
and J is large enough no droplet can appear in the ground state, i.e. the ground state is 
ferromagnetic. In fact, let us assume by contradiction that the ground state is not ai = ~l and 
that there is at least one droplet d in the ground state. Then the energy needed to reverse all the 
spins in 6 must be positive: —"^J^y^YiS) \^\ + ^dip(<5) > 0. Proceeding as above in the proof of 

(4.30) and using that p> d+l,we get Edtp{5) < 'E-yer{S) ItI E|„|>i R? < <^onst. E7er(5) ItI; 
hence for J big enough we get a contradiction and this proves that the ground state is 
ferromagnetic. 

The contour method implemented in the proof of Theorem 3 also allows one to get some 
informations about the geometry of contours and droplets in the ground state for d < p < d+1. 
In fact if 'D(A) is the set of droplets in the ground state and we consider a configuration a' on 
{±1}^ with A' s.t. X>(A') = X>(A) \ {S}, for some 5, then it must be HN{q!) - Eo{N) > 0. 
With the definitions introduced above we have 

Hn{s!) - Eo{N) < -2 J J2 1^1 + ^d^p{S) < 

( pd|Sd|2P-^ . ^^ V- , , 2rf|Srf| |(5| 

<A-J^ '^\_^ +^pffl) E H + 7^^- (4.32) 

7er(i5) 

Imposing that the r.h.s. of (4.32) is positive and choosing ^p(^) + ^ _^ J = —1, we find 

ET,Gr(5) I7I < consi.|(5|e"'^/l^''-il if p = d+1 and E^eri^) \l\ < const. \b\.r^S^ iid<p< d+l. 
US does not wind up Ajv, using the isoperimetric inequality \5\ < (E'verfa) \l\/'^^)'^ ^^ ^^^^ find 

that in the ground state E7er((5) 1^1 — const. e'-'^^^''-^^ iip — d+1 or J^-yeriS) ItI — const.J '^+^-p 
\ld<p<d+l. 

Finally we mention that exploiting the same energy arguments above, one can prove that 
the ground state of Hamiltonian (1.1) with An an A^"^^^ x w cylinder with periodic boundary 
conditions is the ferromagnetic ground state, provided w < K'^e"'''^''-^', if p = d + 1, ov w < 

K' J'i+^-p ii d < p < d + 1, for a suitable constant K' . This result strongly suggests that for 
d < p < d+1 the droplets in the ground state are quasi-((i— 1) dimensional structures of width 

(9(g-^/|S£i-il) or 0(J<i+i-p), reminiscent of the conjectured stripes. If we could prove that the 
droplets in the ground state must necessarily be stripes, then we could exploit the methods in 
Section. Ill to prove that the ground state is indeed realized by a periodic striped configuration. 



V. CONCLUDING REMARKS 
A. One dimension. 

Ground state. In ID the characterization of the infinite system ground state is complete: 
the ground state is periodic and there is a sequence of transitions from an antiferromagnetic 
state of period two to 2/i-periodic states of blocks of sizes h, h > 1, with alternating sign, the 
size h changing (increasing discontinuously) when the ferromagnetic interaction strength J is 
increased. 

Our proof is based on a reflection positivity argument which relies heavily on the details of 
the model Hamiltonian, e.g. on the fact that the long-range repulsion is reflection positive. 
It is an interesting open problem to establish more general conditions the long range repulsive 
interaction should satisfy in order to guarantee existence of a periodic modulated ground state. 
Note that the problem of determining the ground state of a spin system with positive and 
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convex potential was solved by Hubbard and by Pokrovsky and Uimin (even in presence of a 
magnetic field) [23, 24]: this means that for J = the assumption of convexity of the potential 
is enough to determine the ground state of (1.1). The proof in [23, 24] was generalized to an 
"almost" convex case by Jedrzejewski and Miekisz [25, 26]. However this proof requires a small 
ferromagnetic interaction 1 < J < 1 + 2^^. It is an open problem to generalize our analysis (in 
the presence of a large ferromagnetic coupling) to the case of a long range convex interaction. 

It would also be interesting to establish what would be the effect of adding a magnetic field to 
the model Hamiltonian. It is expected that in the presence of a large short range ferromagnetic 
coupling and of a positive magnetic field B of increasing strength the ground state is still 
periodic with the + blocks larger than the — 's, at least if |i3| is not too large [18]. However for 
large magnetic field it could be the case that small variations of the magnetic field could induce 
an infinite sequence of transitions between periodic states characterized by different rational 
values of the magnetization (Devil's staircase): this is in fact what happens for J = as a 
function of the magnetic field [33] . According to [33] the antiferromagnetic state would remain 
the ground state for |5| < Be ior J — 0. For larger values of B there is a sequence of transitions 
from the antiferromagnetic state of period two to (complicated) periodic states with all possible 
rational values of the magnetization. 

Positive temperature. Another interesting open problem in one dimension consists in estab- 
lishing properties of the infinite volume Gibbs measures at positive temperatures. Note that 
for a ferromagnetic long range interaction with decay l/r^, 1 < p < 2, it is known [34, 35] that 
there is a phase transition for the inverse temperature f3 large enough. In the case of model 
(1.1), on the contrary, it is natural to expect a unique Gibbs measure for any finite /3, even for 
1 < p < 2. It is in fact known that if J = in (1.1), then for d = 1 and p > 1 there is a unique 
limit Gibbs state at all values of the inverse temperature /3 [27, 36, 37]. The proof in [27] does 
not extend to the case J > 0. It is known however that for J > all Gibbs states are transla- 
tion invariant [36]. Similarly, using the argument in [38], one can show that the Gibbs states 
obtained as limits of finite volume Gibbs measures with boundary conditions t'^ct (ft,*(J)), 
where g_ {h* (J)) is a periodic ground state configuration with blocks of size h*{J) and t the 
translation operator, are all equivalent among each other, for any k — 1, . . . , h*{J). If, on the 
contrary, < p < 1 and, say, J = 0, it has been conjectured [39] that model (1.1), modified 
as explained after (1.1), admits at least two different Gibbs states obtained as limits of finite 

volume Gibbs measures with boundary conditions ■ • • H \ — • • • and • • • — I \- ■ ■ ■ (which 

are presumably well defined because the Hamiltonian is thermodynamically stable, see Remark 
(1) after Theorem 2). 



B. Higher dimensions. 

Ground state. In two or more dimensions we have shown that the specific ground state energy 
agrees asymptotically for large J with the best variational ground state known so far, which is 
a periodic striped configuration [14]. 

As already mentioned in the introduction, the ground states are believed, on the basis of 
variational computations and MonteCarlo simulations, to consist of a sequence of "stripes" of 
size h (growing when the ferromagnetic interaction strength J is increased) and alternating 
signs. 

If J is sufficiently large, the addition of a magnetic field is believed to lead to a thickening 
of the stripes with spins parallel to the external field. At very large fields a transition from a 
striped to a "bubble" phase is expected [18]: the bubble phase consisting of large "cilindrical" 
droplets of spins parallel to the external field arranged in a periodic fashion and surrounded by 
a "sea" of spins of opposite sign. 

Positive temperature. The non trivial structure of the set of ground states is believed to 
have a counterpart at positive temperature: in particular in the absence of magnetic field it is 
believed that at low temperature there are different pure Gibbs states describing striped states 
with two possible orientations (horizontal and vertical). The striped states are expected to 
"melt" at a positive critical temperature [6] with the "stripe melting" described in terms of an 
effective Landau-Ginzburg free energy functional [2]. 

It goes without saying that it would be of great interest to substantiate these pictures by 
rigorous proofs. Since most standard methods (cluster expansion, correlation inequalities) seem 
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to fail in giving any useful information for systems of spins with long range antiferromagnetic 
interactions, new ideas are needed to understand some of the aforementioned problems. 



APPENDIX A: A PRIORI BOUNDS ON THE SIZE OF THE BLOCKS 

In this appendix we prove (3.13). Given a block _B^. of size hi, let us assume that hi = 
{2K + l)h, for two integers K and h to be chosen conveniently (it will be clear from the proof 
below that this assumption is not restrictive). To be definite let us assume that the block Bh^ 
of size hi we are considering is a + block. Since hi = (2K + l)h, we can think of B^^ as a 
sequence of three contiguous blocks of sizes Kh, h and Kh respectively, to be called Bj^, B^, Bf^. 
Let us now compute the energy needed to flip the block B^ and let us call it AE'{h). If we are 
in the ground state, then of course AE'{h) > 0. Also, we have AE'{h) < 4 J — 2Ei{h), where 
Ei{h) is the dipole interaction energy between a + block B^ of size h and a configuration of 
spins in Z \ Bf^ such that the spins in B\, Bf^ are all + and the spins in Z \ Bf^. are all — . Ei (h) 
is readily computed as 

h h+Kh 
Ei{h) =2^^ ^ '^ 



Z^l Z^ Z^ ^\{j-i)P 



i=l j=h+l j>h+Kh+ 
h h+Kh 

E[E - E 

i=l j=h+l j>h+Kh+l 



_^ckP-ie-"0-») 







r(p) 



so that, using the condition AE'{h) > 0, we find that in the ground state 

W-^ ""' \i ''-M2 (l - e""')(l - 2e-"^') < J (A.2) 

r(p) (i-e °'y 

Using now the fact that ae^"/^ < 1 — e^" < a for a > 0, the l.h.s. of (A.2) can be bounded 
below by: 

t4 r ^^^""(1 - ^""')(l - 2e-"(^''-^)) (A.3) 

T{p) Jq a a^ p 

Now, if p = 2 (A.3) can be computed to give \og{h + 1)/(1 + l/K)'^. Then, choosing K = \ 
and plugging these bounds back into (A.2), we find that, if p = 2, log/i < J + 2 log 2 and 
h, = 5h< 12e"'. 

If 1 < p < 2, (A.3) can be further bounded below as 

1 f^^ rlr\i 1 



-|e-"'*(l - e-"'*)2 - e""^'*(l - e""'*)! 



T{p) Jq a o? 



> 



Choosing K in such a way that B' jK^ ^ < A' /2 and plugging the bounds into (A.2), we find 
that A'h^-P < 2J, so that h, = {2K + l)h < {2K + 1){2J/A'Y/'^^-p\ 



APPENDIX B: THE INTEGRAL REPRESENTATION 

In this Appendix we want to give the details of the computation leading to (3.7). First of 
all let us note that the first two terms in the r.h.s. of (3.7) come from the short range FM 
interaction. In order to show that the computation of the dipole interaction energy leads to 
the remaining terms appearing in the r.h.s. of (3.7), let us first note that the infinite volume 
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dipolc interaction energy £^12 between two blocks of signs £1 and £2, ^a — il, of sizes hi and 
/i2 and separated by a string of d spins, can be computed as follows: 

1 / aa p-i a{j-i+hi+d) _ 






Now, in order to compute the finite volume dipole interaction energy of two different blocks hi 
and hj, j > i, we can simply apply the definition and (B.l) to find: 

h^ hj 

^1-^1 2^ 2^ 2^ Itti-Z-U/i-J Uh. , ^O-nNlP ^ \ ■ ) 



-^ 3 /L^ z^ ^L^ \m-l + hi^ h /i,-i + 2nN\P 



nGZ / 

/•oo 



-2anN 



•'" ^^' ^ ' i<k<j n>a 



M 



which leads to the term in the third line of (3.7) (we used that SiSj ~ (—1)"' * and 'Ylii=\ ^i — 
2N). Similarly, the finite volume dipole self-interaction of a block hi can be computed as 

v' ^'-fe I 2V V ^ - r "^^ a^-^ ^'° fB3) 



■ k(l - e"") - (1 - e-"''') + 2e"'**(l - e""'**)^ 



1-e 



-2aN 



Summing (B.3) over i — 1, . . . , M, we get the last term in the first line of (3.7) and the term in 
the second line. 



APPENDIX C: COMPUTATION OF THE ID ENERGY GAP 

In this Appendix we want to compute the specific energy e of the [Ah* + 2)-periodic config- 
uration (ft,*, /i* + 1, ft* + 1, /i*, . . .) in the case in which ft* and ft* + 1 are both minimizers of 
e(ft). In particular we will show that e — cq = 5 > Q. Using the general expression (3.7), after 
some algebra we find that 

. . 4J 2 r, , j2(l-e-"''*)(l-e-"(''*+i)) 



2ft* + 1 2ft* + l7o 1 l_e-a(4h*+2) 

^-2ah' M _ g-a(h*+l)\2 ^ ^-2a{h' + 1) ^^ _ ^-ah'\2 ^ 
I _ g-a(4h*+2) f 



(C.l) 



where i^{a) = {T{p)) ^a^ ^e "(1 — e ") ^. We want to prove that this expression is strictly 
larger than gq. Note that, since both ft* and ft* + 1 are minimizers, we have 

— / da z^ a tanh = - / daz^a tanh C2 

ft* ft* io 2 ft* + 1 ft* + 1 io 2 ^ ' 



implying that 



J= / dau{a) (ft* + 1) tanh ^ - ft* tanh "(^* + ^^ (C.3) 



and 



'0 

Using (C.3) and (C.4) we find 
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~ / , . / ah* a(h* + 1)\ , , 

eo = 2 / da ly (a) [tanh-- tanh ^ -) (C.4) 



e — Cq — 2 i da v{a) 







, ah* , a{h* + 1) ,„ , 

tanh — - + tanh — -— '- - (C.5) 



2(1 - e^"'**)(l - e""('''+l)) + e^^"'''(l — g-O'ih'+'i-)^ ^ g-2a(h* + l)^ _ p-Qh*\2 

1 _ g-Q(4/i*+2) 

A bit more of algebra shows that (C.5) can be rewritten as 

f-OO 

6-60 = 2/ rfai.(a)(6-"'** -e-"('''+i))2(l-e-"''')(l-6-"(''*+i)) (C.6) 

and this concludes the proof. 
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